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Abstract 

This report describes a dynamic gas temperature measurement system. It has frequency 
response to 1000 Hz, and can be used to measure temperatures in hot, high pressure, high 
velocity flows. A personal computer is used for collecting and processing data, which results 
in a much shorter wait for results than previously. The data collection process and the user 
interface are described in detail. The changes made in transporting the software from a 
mainframe to a personal computer are described in appendices, as is the overall theory of 
operation. 


Introduction 

The dynamic gas temperature measurement system (refs 1-4) consists of a probe with two 
thermocouples of different diameters (fig. 1) and special software. The signals from the 
thermocouples are fed through differential amplifiers and into the data collection/digitizing 
board in a PC. There are 3 signal channels. The first two are AC coupled and are the 
fluctuating part of the signal from the two thermocouples. The third channel is DC coupled 
and is the mean temperature signal from the large thermocouple. By comparing the 
amplitudes of the two signals at low frequencies, the software is able to provide the dynamic 
gas temperature, frequency compensated to 1 kHz. A brief description is given here; the full 
details are in appendix A. 

The software begins by calculating the frequency response function for each thermocouple 
using the one dimensional model shown in figure 2. Because the thermocouples are bead- 
less, they can be treated as cylinders in cross-flow, so that the temperature distribution in the 
wire obeys the one dimensional heat transfer equation: 


dT _ d 2 T 
Tt ~ 01 dx 2 




( 1 ) 
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where T = T{x,t) = wire temperature at location x & time f, K 
a = wire thermal diffusivity, m 2 /s 
p - wire density, kg/m 3 
c = wire specific heat, J/kg/K 
D = wire diameter, m 
T t (t ) = instantaneous gas temperature, K 
h = convective heat transfer coefficient, W/m 2 /K 


The first term on the right is required because the low length to diameter ratios of the 
thermocouples require that heat conduction down the wire be accounted for in calculating the 
wire temperature. 

The frequency response depends on the heat transfer coefficient h as well as on D, the wire 
diameter. The value of D is known, and based upon user input and the measured mean gas 
temperature, the program makes an initial estimate of the value of h. Actually, a related 
parameter, T, is used because it is independent of diameter and has the same value for both 
thermocouples. It then goes on to calculate the ratio of the frequency response functions for 
the two thermocouples for values of T ranging from 0.2 of the estimated value to 1.8 times 
the estimated value in steps of 0.1. Each of the calculated ratios, together with the 
corresponding value of T is then stored, for later comparison with the measured ratio. 
Determining measured frequency response ratio is the next step to be performed. 

The temperature data consists of the time varying part of the signal from each thermocouple, 
and the mean value from the large thermocouple. Typically in a run, about 15 seconds of 
data are digitized and recorded. This time record is then broken into blocks of 1/2 second 
each, and each of the blocks is Fourier transformed. Next, the auto and cross spectral 
density functions are calculated from sums of the transformed data blocks, and finally, the in 
situ frequency response ratio is obtained from the ratio of the cross spectral density to(8n)of 
the auto spectral densities. The program then compares this ratio to the set of ratios q'/v^ 
previously calculated, in order to determine the corresponding value of T, interpolating if 
necessary. The frequency response then calculated using this value of T in the mathematical 
model is taken to be the true frequency response of the thermocouple under the prevailing 
flow conditions and is used to frequency compensate the thermocouple signals. 

As originally implemented, data reduction for dynamic gas temperature measurements 
required the services of a central computing facility and a mainframe computer. The time 
interval between the collection of data and the availability of results could be as long as a 
week. 

Due to the availability of more powerful personal computers, and to hardware which permits 
rapid collection and storage of data on the hard disc, it is now possible to collect and reduce 
data using a personal computer located at the test facility. The results of a measurement are 
now available in minutes instead of days, and the process is more user friendly. 
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In transferring the Dynamic Gas Temperature Measurement System to a personal computer, 
several changes were made to the software. The first was to replace the numerical technique 
used to calculate the frequency response function with a closed form solution. This one 
change resulted in a reduction of execution time from about 7 minutes of CPU time to about 
90 seconds. The derivation of the formula is in Appendix B. The second change was to the 
fast Fourier transform routine, which resulted in a further reduction in execution time to 23 
seconds. There are at least 60 transforms to be performed, 30 for each thermocouple, and 
each transform is 2048 points. By making use of the fact that the data consists of real 
numbers rather than complex numbers, it was possible to cut in half the execution time for 
performing a transform. The method works by treating the 2048 real data points as 1024 
complex data points, performing the transform, and unscrambling the results. The method is 
described more fully in Appendix C. 

In addition, major changes were made to the way the program accepts user input, and the 
way data is acquired from the thermocouples. These changes are described in the following 
section. 

Software 


The data acquisition specification for the Dynamic Gas Temperature Measurement System 
requires the acquisition of three channels. These channels are the AC part of the signal from 
the small diameter thermocouple, and the AC and DC parts of the large diameter 
thermocouple, called respectively: small thermocouple-AC, large thermocouple-AC, large 
thermocouple-DC. These three signals must be simultaneously sampled at intervals of 
approximately 244.1 micro-seconds for a duration of up to 60. seconds. 

Performing this task is a Data Translation DT2829 Acquisition Card used in conjunction with 
a Compaq 486/25 Personal Computer. Some important features of the DT2829 are the eight 
single-ended analog inputs with simultaneous sample-and-hold capability, the 16-bit A/D 
converter operating at a throughput rate of 30Khz, the external clock and trigger options and 
the dual DMA transfer capability which allows continuous acquisition to disk. 

The program written to use the Data Translation Acquisition Card is graphically based and 
menu driven. The user begins the program by typing Acquire. Once inside the operating 
environment, the user selects OPEN from the list of options in the main window menu 
named FILE (figure 3). After opening a file for the data, the user can start filling in 
acquisition parameters such as the clock source, the scan time interval and the channel 
selections (figure 4). 

For each of the eight channels there is a BNC connector icon. There also is a connection 
point for each of the eight locations in the scan list. The DT2829 works its way down the 
scan list reading the channels assigned to each scan list location. Any channel can be in any 
scan list location or in several scan list locations, but the first channel is usually defaulted to 
the first scan location because the hardware samples-and-holds the current values in the other 
seven channels upon each acquisition of the first channel. In order to connect a particular 
channel to a particular scan location, the user simply draws a wire starting from the desired 
scan location to the desired BNC channel icon. Afterwards, the user is prompted for an 
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external gain and a descriptive comment for that particular channel. If the user desires to 
change the gain or the comment fields, the cursor is clicked in either field and a prompt for 
new information mil appear. The scan list must be filled in ascending order to prevent the 
user from creating voids where a scan location has no associated channel. In order to 
remove a channel from the list, the user clicks the cursor twice at the associated scan 
location. 

The requested acquisition time duration and requested time interval between acquisitions can 
be changed by once again clicking the cursor in these fields and responding to the prompts. 
After entering the requested time duration and time interval, the actual values for these 
quantities are displayed in an adjacent field. If an illegal value is entered, the user will be 
notified. The clock mode (internal or external) is selected by the clicking the cursor in the 
desired field. Selecting the external clock mode will cause the program to prompt the user 
for the external clock speed. Similar to the clock mode, the trigger mode (internal, external 
and external scan) is selected by clicking the cursor in the desired field. The user is 
prompted for period information when the external-scan trigger mode is selected. 

Once the data storage file is chosen and all of the parameters are setup, the user can acquire 
data by selecting ACQUIRE from the list of options in the main window menu named 
OPTIONS. After completing the acquisition, the data can be displayed by selecting DISPLAY 
from the same main window menu named OPTIONS. 

All recorded channels can be viewed individually and scaled along the time axis by using the 
appropriate display window menu options (figure 5). When satisfied that the data is suitable, 
the user exits the display window. Afterward, the user can either close the data file and 
request that it be saved or not saved, quit without saving the data file or quit and save the 
data file. Any of these choices can be selected from the main window menu named file. All 
saved files will have the form ’Filename, acq’ and can not be re-displayed with the program 
Acquire because the data file will be erased upon re-entry into the program. 

The data files created by the program Acquire are read by another program called Convert 
which allows the user to reconfigure the data into a text file, an IEEE real file or a data file 
compatible with the Dynamic Gas Temperature Program. The user simply types Convert to 
enter the program and then selects an acquisition file using OPEN from the main window 
menu named FILE (figure 6). Then the conversion type (ASCII, BINARY or ANALYZE) is 
selected from the main window menu named OPTIONS. Individual channels are selected for 
conversion by clicking the cursor on the desired conversion output field and then clicking the 
cursor at the desired input channel field (figure 7). Once all of the desired channels have 
been selected for conversion, select CONVERT from the main window menu named 
OPTIONS. The user will be asked if the output data is to be divided by the external gain for 
an actual reading or kept in its original form. After converting the data, the program is 
terminated using the QUIT option from the main window menu named OPTIONS. The 
resulting output file will be in one of three forms (Filename.txt, Filename.bin or 
Filename.dat) determined by the conversion type selected by the user. If the data file for the 
Dynamic Gas Temperature Program was created ( Filename.dat ), the menu program Thermo 
can be used to access its information. 
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The original program was written in FORTRAN, a language designed with rudimentary 
Input/Output capabilities. Because of this, the I/O section of the data reduction software is 
cumbersome to use. 

Rather than rewrite the I/O section of the FORTRAN code to allow for the improved I/O 
capabilities, we decided to develop a program which would become the input interface to the 
FORTRAN program from the user. This code, written in Borland Turbo Pascal Version 5.5, 
is called Thermo, and is described here. 

After acquiring and converting the data, and assigning a data set name, the user calls 
’Thermo’ to set up the input file for the data reduction software. The first concern of the 
code is to find and list all the data sets found on the disk. While it is doing this, the 
program displays the message shown in fig 8. After the files are sorted, the Main Menu is 
displayed. With the pressing of the ’Enter’ key, the Main Menu screen appears as in fig. 9. 
The Directory is shown in the bottom half of the screen, in alphabetical order. Those files 
marked with an asterisk are files for which the data reduction software has already been run. 
After choosing a data file, the user must decide to either run the FORTRAN program 
{EXECUTE)-, edit the input {MODIFY), quit the program {QUIT) or ask for help {HELP). 

If the user requests help in the Main Menu, the screen displays the contents of figure 10. 

This screen displays the available functions for the program, along with the keystrokes which 
activate them. 

If the user edits the input file, a full-screen data editor is invoked. Pressing ’Help’ (F3) 
from the data editor displays fig. 11. Fig 11 shows the data to be modified, and each entry 
is individually modifiable. After accepting any changes, the data editor generally takes the 
user through a series of questions to finish the data entry process. These questions are 
variable, based on the entries in the first part of the data editor. If the user wants to save the 
modified input file, the program will require verification of the file name. This is simply a 
safeguard to ensure that the proper data set was edited. The stored input file is functionally 
equivalent to the file described in reference 2. This is the file used by the processing 
software to determine the data manipulation and display required by the user. The seven 
possible display outputs, corresponding to the 7 example display types in ref. 2 are shown in 
Appendix D. 
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Figure 1: Dual Wire Thermocouple Probe 
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figure 2: One Dimensional Model of Thermocouple. The model is symmetric about the 
junction. 



Figure 3: Data Acquisition Window — Menu Choices 
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figure 6: Data Conversion Window - Menu Choices 



Figure 7: Data Conversion Window - Window Format 
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SORTING DATA FILES 


Figure 8: Initial Screen; Waiting for Sort to Finish 


Dynamic Gas Temperature 
Measurement System 
Interface 

NASA Lewis Research Center 
Drive X: No File Selected 

MAIN MENU 


DIRECTORY 

EXECUTE 

MODIFY 

QUIT 

HELP 

Directory of X: 




FILE1 

FILE2 

FILE3 

FILE4 

FILES 

FILE6 

FILE7 

FILE8 

FILE9 

FILEO 


Figure 9: Main Menu with Directory Option Chosen 
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Dynamic Temperature Measurement System Available Functions: 


FI 

F2 

F3 

F4 

F5 

F6 

F7 

F8 

F9 

KM 

Quit 

Edit 

Help 


File 


Exec 



Em 


F1- 

F2- 

F3- 

F4- 

F5- 

F6- 

F7- 

F8- 

F9- 

F10- 


Quit - Exit w/o Save/Execute. 
Edit - Edit a file. 

Help - Show this screen. 

File - Get Input File. 

Exec - Exit & Save/Execute. 

FFT - Delete FFT files. 


Esc - Escape - 
A C - Ctrl C - 
Horn - Home - 
End- End - 
PgU - PgUp - 
PgD - PgDn - 
Ins - Insert - 
Del - Delete - 
BkSp-BackSpce- 


Exit w/o Save/Execute. 
Exit w/o Save/Execute, 
beg of entry/ line/screen, 
end of entry/ line/screen. 
Toggle betwen nun/text. 
Toggle betwen nun/text. 
Toggle insert mode 
Delete character. 

Delete prev. character. 


Files list the data files on C:, and load a saved file. 

Execute runs the analysis software using your selected file 

Edit will allow you to edit a loaded file or create a new one 

Quit terminates the program 

, Press any key to continue: 


Figure 10: Main Menu Help Screen 


TC0(1,1) TC0(2,1) TCD(3,1) TCD<4,1) TCD<1,2) TCD<2,2) TCD(3,2> TCD(4 t 2> 

GAS(I) GAS(2) GAS(4) GAS(3) BLS2 >(1) -(2) 

FREQ<1) FREQC2) FREQC3) FREQ(4) PLTFRQ * CHAKL(I) CHANL(2) 

CHAN L (3) CHANL(4) CHANL(5) CHANLC6) CHANLC7) CHANL(8) CHANL(9) U 

TEXT DATA GOES HERE 


KEY: - = IAVDAT, * = 1FLAGSC2), # = PltChc 
See data editor for more information. 

Press any key to continue: 

Figure 11: Help Screen Available in the Editor 
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Appendix A 

Theory of Operation of the Dynamic Gas Temperature 
Measurement System 


For purposes of computation each of the thermocouples is represented by the one 
dimensional model shown in figure 12. The thermocouples are constructed without beads, so 
can be treated as cylinders in crossflow. The temperature distribution in the thermocouple 
then obeys the one dimensional heat transfer equation 


dT 

dt 


a— + —(T(t)-T) 
dx 2 pcD * 


(A.l) 


where 

a = thermal diffusivity of the wire, m 2 /s 

p = wire density, kg/m 3 

c = wire specific heat, J/kg-K 

D = wire diameter, m 

T g (t) « instantaneous gas temperature, K 

h = heat transfer coefficient, W/m 2 -K 

The first four quantities are known, either from handbooks or by direct measurement. The 
gas temperature, T g (t), is what we are trying to measure, and is not known. And the value 
of h, the heat transfer coefficient, since it depends on local flow conditions, is not known 
either. However, although the exact value of h for the particular operating conditions is not 
known, it can be estimated based on the mean temperature from the large thermocouple (ref. 
5). 

Having chosen a starting value for h, the thermocouple is treated as a linear system (ref. 6). 
From linear system theory, the response to a sinusoidal input e 2,rift is H r (f)e 2lrift where H r (f) 
is the frequency response function. The subscript.!! refers to the parameter 

r = hH L < A * 2) 

pc 
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which is used instead of h because it has the same value for both thermocouples. 

The next step, is to use the mathematical model of the thermocouple to generate the ratio of 
frequency response function magnitudes | H r iQ (J) | / 1 H r Jf) | for values of T from 0.2r e to 
l-8r e , in steps of 0. ir e . The quantity I\ is the value of T based on the estimated value of h. 
The subscripts 3 & 10 refer to the diameters of the thermocouple wires. This set of ratios is 
then saved, to be compared with the measured ratio obtained from the data. 

To measure the value of T in the flow, a time history of length t r seconds is recorded from 
each thermocouple. It is assumed that flow conditions are steady during this period. The 
thermocouple data is linearized using a polynomial appropriate to the thermocouple material 
(ref. 7). The thermocouple can now be regarded as a black box producing an output y(r), 
which is directly proportional to the input x(t), the gas temperature. 

The t x seconds of data are separated into rz, records, each of length t seconds, so that nj = t r . 
Typically 15 seconds of data is recorded, which is broken up into 30 blocks of l A second 
each. Breaking up the total time record into small pieces tends to improve the statistics (see 
for instance section 8.6 of ref. 6) 

The Fourier transform of each record is now computed numerically 

n.,00 - «TM'>] (A.3) 

k = l,2,-», n r 


where FFT stands for fast Fourier transform. Note that capital letters are used for frequency 
domain data. 

The auto and cross spectral density functions (ref. 8, sections 3.3.2 
and 3.3.4) are then calculated from the Fourier transforms: 


G »W - ^E I y «WI 2 

n r L *«1 

g mo« - ^E Km„<S) 


(A.4) 
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The reason for calculating these quantities is that they are related to the frequency response 
functions of the thermocouples via the relations 


G„W • 

G,., ' XWiJfiGJA 


The quantity G„ (/) in the above equations is the unknown auto spectial density of the gas 
temperature fluctuations. 


The measured gain rado of the two thermocouples is now found by merely dividing, i.e., 


g 3jo 0 _ h;w k w 

G 3> 3 (f) " \HM 2 


\m\ 


((cos^.3 -<#> 10 ) +zsin(0 3 -4> 10 )) 


(A.6) 


and the desired ratio | H l 0 (f) | / 1 H 3 (f) | is just the square root of the sum of the squares of 
the real and imaginary parts of G 3 AO (f)/G 3 > 3 (f). 


The measured value of the gain ratio is then compared with the set of gain ratios calculated 
earlier. Once a match is found, the in-situ value of T has been found. This value of T is 
used with the mathematical model of the thermocouple to generate the frequency response 
function for the small thermocouple, H 3 (j). Then the Fourier transform of the gas 
temperature fluctuation, X(j), is related to the Fourier transform of the signal from the 
thermocouple, Y 3 (f), by the relation 


m 


y 3 (J) 

m 


(A. 7 ) 


Finally, the unknown gas temperature is just 

r,(f) = x(o = bt-'[x(/)] 


December 14, 1992 


14 


Appendix B 

Frequency Response Function 


This appendix contains a derivation of the analytic expression for the frequency response of a 
supported thermocouple wire in which there is conductive heat transfer down the wire (fig. 
12). With reference to the photograph of the probe shown in figure 1, the smaller wire in 
figure 12 represents the horizontal cross wire, & the larger diameter wire represents the 
vertical support wire, part of which is exposed to the flow and part of which extends into the 
ceramic probe housing. The expression derived in this section replaces the finite difference 
method originally used. 

The temperature in the thermocouple wire, T w (x,t), satisfies the one dimensional heat transfer 
equation 

= (BJ) 

where p = wire density, kg/m 3 

c — wire specific heat, J/kg-K 
k = thermal conductivity, W/m-K 
h = convective heat transfer coefficient, W/m 2 -K 
D = wire diameter, m 
T t = gas temperature, K 

In our case, the gas temperature is assumed to be a function of time only, and to be 
described by a sinusoidal fluctuation superimposed on a mean value, 

T g (t ) = T+T < e u * (B.2) 

where T — mean gas temperature 

T a = amplitude of sinusoidally varying part 

The wire temperature is then also the sum of a time independent and a time dependent part. 

Let 

TJx,t) = TJ$+T{x,t) (B.3) 


and substitute equations (B.2) and (B.3) into equation (B.l): 


dT 

w 

~dt 


dT 

Bt 


a- 


d 2 T„ 

dx 1 



w'here a = thermal diffusivity of the wire, m 2 /s. 
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This equation can be separated into two, one of which is independent of time and one which 
is not: 

c£^.(T-7-J . 0 <B.4> 

■ 

dx 2 V ' 

It is equation B.5 which is of interest here. The steady 
obtained by separating variables. Let 

T(x,t) = X(x)e iul 

and substitute (B.6) into (B.5). After cancelling the common factor and gathering terms, 
the result is 


££ (B.5) 

dt 

state solution of equation (B.5) is 

(B.6) 


. . o) 
1+z — 


X - - — +T 

<*:dx 2 - 


Upon dividing through by T a to normalize, the result is 

-Giu>)f = -1 


(B.7) 


where 



G = l+UL 


(B.8) 


Equation (B.7) is an ordinary differential equation which has the solution 

fix) = /4cosh(<jx) +£sinh(<jx) +-^ 

G 


(B.9) 
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where G is defined in equation (B.8), q = ( Gly ) m , and A and B depend on the boundary 
conditions. 


Thus for the three regions shown in figure 12, 

fj = AfOsh(q/) +B j sinh(q J ic) + JL, y-1,2,3 

'J/ 


except that 1/G 3 = 0. That is because the part of the support wire in region 3 is not exposed 
to the flow, so h 3 = 0. Then o> 3 = 4 h 3 /pcD = 0 and 


i -i 



However, 


«? - 3 - 

0> 3 


f * 

1+/-5L 


a 


o > 3 


The unknown constants A i & B } are determined by the requirement that 

a) the temperature is continuous everywhere 

b) the heat flux is continuous everywhere 

c) the temperature fluctuations at x=l 3 are zero, i.e.,/ (1 3 ) = 0 

d) the material properties are the same on either side of the junction, so the solution is 
symmetric about x=0. 


(This latter simplification of uniform material properties is adequate for the type B 
thermocouples used in high temperature environments, but not for thermocouples such as 
type K used for lower temperatures. In such cases the solution given here must be extended 
to take into account the differences in material properties, as has been done in refs 9 & 10). 

Condition d) states that/ (-*) =f (x), which in tum implies 


so that 


f(0) = (^sinh^) +<j 1 B 1 cosh(< 7 r r)) 


x=0 


= q l B l = 0 
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(B.10) 


/,(*) = i4,cosh(9 r x)+-L 

G i 

Because of the boundary condition c) on/ 3 it is convenient to write / 3 as 

f 3 (x) = A 3 cosh(q 3 (l 3 -x)) +B 3 sinh(^ 3 (/ 3 -*)) 

Then f 3 {l 3 ) = A 3 = 0, and 

f 3 (x) = 2*3 sinh(<7 3 (/ 3 -x)) (B.ll) 

The six unknowns A u B u /1 3 , B 3 have been reduced in number to four, which will now be 
found by applying the requirements a) & b) at x=/ 1 & at x=l 2 . Continuity of temperature 
and heat flux at x=l u 

fA) = / 2 ft) 

and at l 2 

m -m 

-/A) 

give four equations in four unknowns: 

/l 1 cosh(q 1 / 1 )+J r = A 2 +-^ r (B.12) 

G l G 2 

^j^sinh^j/j) = R q^ 2 (B.13) 

^cosh(^ 2 (/ 2 -/,)) + £ 2 sinh(tf 2 (/ 2 -/,)) = ^ 3 sinh(^ 3 (/ 3 -/ 2 )) (B.14) 

v4 2 ^ 2 sinh(^ 2 (/ 2 -/j)) + B 2 q 2 cosh(q 2 (l 2 -l l )) = -A 3 q 3 cosh {q 3 (l 3 -lj) <®- 15 ) 

where R = D 2 /D ? andf 2 {x) has been written in the form 

f 2 (x) = ^ 2 cosh(^ 2 Cx-/ 1 ))+B 2 sinh(g 2 (x-/ 1 ))+-L 

°2 

The desired quantity is the temperature at the thermocouple junction, which from equation 
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(B.10) is/i(0)= A x + 1/Gj. The quantity G, is defined in equation (B.8), and A x is extracted 
from the above equations (B. 12) to (B.15). The result is 


A, = Rq r 


g 2 "c?, 


q 2 sinh(q 3 (l 3 -lj)swh(q 2 (l 2 -l i )) 


Rq 2 


^cosh( ?3 (/ 3 -/ 2 )) 


J_ _1_ 


cosh 


(B.16) 


where 


A - ^ 2 sinh(^ 3 (/ 3 -Z 2 ))^?g 2 cosh(^,/ 1 )sinh(^ 2 (/ 2 -/ 1 ))+^ 1 sinh(^ 1 / 1 )cosh(^ 2 (/ 2 -/ 1 )j) ^ 1?) 

+ 5 3 cosh(^ 3 (/ 3 -/ 2 ))(^ 2 cosh(^ 1 / 1 )cosh(g 2 (/ 2 -/ 1 ))+^ 1 smh(^ 1 / 1 )sinh(^ 2 (/ 2 -f 1 ))) 

Equations (B.16) and (B.17) form the basis for the changes made in subroutines TRFP and 
TRFEM. They replace the finite difference equations formerly used. 
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Appendix C 
Modified FFT Algorithm 

Typically, there are 30 blocks of data taken for each run, with 2048 data points per block, 
and each of these blocks of data is Fourier transformed. Since there are two thermocouples, 
there are 60 transforms for each data run. Because of the large number of transforms to be 
calculated, an increase in computational efficiency could significantly reduce execution time. 


Let a typical data point in one of these blocks of data be written x(k), where £=0,1,2,--,2AM 
is a free index denoting location within the data block. The size of the data block is 2 N, 
which is a power of 2 (ref. 11). In our case 2 N = 2048. 

Fast Fourier Transform programs are designed to calculate the transform of an array of 
complex numbers. Each of the x(k) must be converted to complex form by writing it as a 
complex number x c (k), whose imaginary part is zero: 

x c (k) = [r(£),o] = x(k)+iO (C.l) 

The transform of the array x c (k) then consists of the array 

m 

n 

and 4 N numbers are calculated, the real and imaginary parts of each of the X(ri). Since we 
started with an array of size 2 N, half the 4 N numbers must be superfluous. This is, in fact, 
the case: 


2N-1 


M WT-J 

J-V x(k)e 

2 NU 


2*ibi 

~TR~ 


= 0, 1, -, 2N-\ 


(C.2) 


First, X(0) is real, which follows directly from the definition (C.2); i.e., with n— 0, the 
expression is just a sum of real numbers, which must be real. Further, since the x c (k) are 
real, the transformed array satisfies the relation 

X(N+m) = X'(N-m), m = 0, 1, ~, AM < c 


where the * denotes complex conjugation. Setting m = 0 in equation (C.3) shows that X(N) 
is real, and for m 5 * 0, die second half the array X(ri) is just the complex conjugate of the 
first half. Thus the Fourier transform X(ri), n=0,l, 2AM, of the real array x(k), £=0,1, 
•••, 2AM, has only the 2 N independent quantities X r (0), A^(l), -, X T (N), X ; (l), X { (2), -,Xi(N- 
1), where r & i denote the real and imaginary parts of X(n). 

Thus padding each real number with a zero as in equation C. 1 is not a very efficient 
procedure for calculating the discrete Fourier transform. Fortunately there is a better way, 
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(refs 11 & 12). By assembling the real array into pairs, and treating each pair as a single 
complex number, one calculates an N point transform instead of a 2N point transform, 
cutting the computing time required roughly in half. The details of how this is done follows. 


As above, given the real array x(k) of size 2 N, the discrete Fourier transform is the array 

2N-1 2 xibi 


.f *4T 1 

X(n) = n = 0, 1, 2N-1 

IN w 


(C.4) 


Grouping the sum into the even and odd numbered parts, 


m 


1 / 2«i2n 

■ 4; k0) + ^2)e‘®' + - 

2N u 


x(2N-2)e 


2*i(2N-2)n\ 

277 


I 2 Tin 

+ yc(l)e r7jr + - +, 


k(2N-l)e 


2nQN-l)n 

w 


)] 


IN 


N- 1 2*ii ffl*» N- 1 2*i(2/+l>i 

E T +E 

j’O j-o 


- AM 2*f» *oi - AM 

iE *&> T + « ""tbE *&> 

/V 7-0 y.o 


2x971 


(C.5) 


where 


h(J) = *(2/) 

= *(2 M) - 

J = 0, 1, - , AT-1 


(C.6) 


The two sums which appear in equation (C.5) are the Fourier transforms of A and g, 

AM 2*9» 


»(») = iE h oy 


7-0 

tf-i 2 x 9B 


N‘ 

m = ^E *(/> 

iVy.0 

n = 0, 1, - , tf-1 


(C.7) 


so 

(Note that the range of the free index n is different in equations C.7 and C.8. This doesn’t 
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xw - 3 


■win 

,-7Ti 


(C.8) 


[H(n)+e "GW J 
n = 0, 1, - , 2N-1 

matter because the functions H and G start to repeat after n exceeds AM). 

The functions H and G are not computed separately, since that would offer no savings in 
computer time. Instead, define the complex array 

y(j) = m+m (C - 9) 


and calculate the single N- point transform 


AM 


2*yn 


m - 

N£o 


2 t fn 


Y(n) = H(n)+iG(n), n = 0,1, - , N-l 


(C.10) 


The trick now is to express H(ri) & G{ri) in terms of Y(ri). First, we need the result that, 
since g(f) and h(j) are real, 

G r {N-n) = G r (n ) (C.I1) 

Gfff-n) = -Gfii) 


and similarly for H. This result is the same as G{N-n) = Cf (n), which follows from 
equation (C.7). Then, since 

Y(n) = H(n)+iG(ri) = p/n) + /Af,.(n)] + /[G r (/i)+/G,.(n)] 


Yfii ) = H r (n)-Gfn) 

(C.12) 

Yfn) = HfiO+GJW 

(C.13) 
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Then, from equations (C.12) and (C.ll), 

Y r (n) = HJrt-Gfn) 
Y r (N-n) = HJrt+Gfr) 



H r (ri) = i(Y r (n)+Y r (N-n)) 

G,(n) ■ -j(W-W-»)) 

(C.14) 

Similarly, equation (C.13) gives 




(C. 15) 

so that 

ff(B) = H,(n)+iH t (n) 

(C.16) 

and 

G(n) = G r (n)*iGfn) 

- ^m*TXN-n))-i(r r (n)-r r (N-ri)) 

(C.17) 

Equations (C.16) and (C.17) can now be substituted into equation (C.8) to give the elements 
of the desired transform X(n) in terms of the downsized transform Y(n). 

The process can also be reversed, so that the original-array *(£) can be recovered from the 
array X(ri) via an efficient inverse transform. After some rearranging, H(n) and G(ri) can be 
written 


H(n)^m+Y’(N-n)) 

(C.18) 


<Xn)~!fm-Y‘lN-n . )) 

(C.19) 


and equation (C.8) becomes 
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(C.20) 


-iwn 

X(n) = 1 to + Y * (N-nj)-le "^to - Y * (N-ri)) 


Then 


X(N-n) = i(Y(N-n)+Y'(n))+le^(Y(N-n)-Y-(nj) 


and 


X*(N-n ) = l(7(/i) + Y * (N-nj) +U to - Y * (N-nj) 


(C.21) 


Adding and subtracting (C.20) and (C.21), 

X(n)+X-(N-n) = i(Y(n)+Y'(N-nj) 

iim 

X(n)-X'(N-n) = ie T (F(n)-F-(AT-n)) 


from which 


ixn 

Yin) = fen) +X * (N-ri))+ie '^(nyx * (A-n)) 


(C.22) 


The inverse iV point transform is then performed on T(n) to give the array y(f), from which 
the array x if) is available via equations (C.9) and (C.6). Equations (C.20) and (C.22) are the 
basis for the subroutine REALFT, which performs all of the Fourier transforms in the 
dynamic gas temperature measurement software. 
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Appendix D 

Plotting Choices available from User Interface - 


These plotting choices correspond to the example plotting outputs found in ref 2. They are 
shown as figures D. 1 through D.7. 
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Figure D.2: Compensated Data: Average Power Spectral Density (K 2 /Hz), 
Instantaneous Temperature (K), and Intantaneous Power Spectral Density (K 2 /Hz) 
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Figure D.3: Compensated Data: Average Log Power Spectral Density (dB), 
Instantaneous Temperature (K) 
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Figure D.4: Compensated Data: Averaged Linear Power Spectral Density (KA/Hz), 
Instantaneous Temperatue (K) 
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RVERRGED FREQUENCY DQMRIN DRTR 


SMRLL T/C 


COMPENSATED DRTR 


BURHER CENTER LINE, U IN. DONNSTREAH 
PROBE « I. REPAIRED <J/BB DRTR RECORDED 


STARTING REC NUMBER 1. 
RECORDS IN AVERAGE 3. 


20. 01 Ub.02 60.03 Bb.05 

FREQUENCY (HZ) *10* 


Figure D.5: Compensated Data: AVeraged Narrowband Power Spectral Density (K 
rms) 
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Figure D.6: Uncompensated Data: Averaged Narrowband Power Spectral Density (K 
rms), Instantaneous Temperature (K), Instantaneous Narrowband Power Spectral 
Density (K rms) 
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COMPOSIT INSTANTANEOUS TIME WAVEFORM 

SMALL T/C COMPENSATED ORTA 

INSTRNTANEOUS DATA, RECORD NUMBER 1 . 

Don.P' CEMTER UNE > i in* downstream 

PROBE « 1. REPAIRED U/BB DflTft REP 



B7.987 K RMS 
1202.123 If MEAN 


0.20 0.30 
time [SEC) 


Figure D.7: Compensated Data: Instantaneous Temperature (K) 
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